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Abstract 

We investigate a particle physics model in a six-dimensional spacetime, where 
two extra dimensions form a torus. Particles with Standard Model charges are con- 
fined by interactions with a scalar field to four four-dimensional branes, two vortices 
accommodating ordinary type fermions and two antivortices accommodating mir- 
ror fermions. We investigate the phenomenological implications of this multibrane 
structure by confronting the model with neutrino physics data. 
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1 Introduction 



The idea that the physical spacetime might have more than three spatial dimensions dates 
back to the beginning of the 20th century and to the attempts of Nordstrom, Kaluza and 
Klein to unify gravity and electromagnetism pQ . These early ideas have been resurrected 
in the present-day string theories which unify gauge theories and gravity in a ten or eleven 
dimensional spacetime. Extra dimensions were originally thought to be compactified to 
extremely short distances. However, according to some new ideas developed within the 
string theories some of them could be relatively large and could have phenomenological 
consequences testable in physics experiments. 

It has been realized that the presence of large extra compact spatial dimensions could 
make the fundamental Planck scale M* , that is the Planck scale associated with the 
higher-dimensional spacetime, much smaller than it is in the ordinary four- dimensional 
world 11] • This would offer an attractive solution to the hierarchy problem, i.e. the dis- 
parity of the fundamental energy scales of weak interactions and gravity. The weakness of 
the gravity as seen in our visible four-dimensional spacetime, or in other words the large- 
ness of the effective four-dimensional Planck scale Mpi ~ 10 19 GeV, could be understood 
to follow from the spreading of the gravity force lines into the extra dimensions. Gravita- 
tional interactions have not yet been tested at short distances, and at present experiments 
still allow some of these extra dimensions to be as large as tenths of a millimeter [HUH]- 
In order to have the Planck scale M* at the TeV range, as the solution of the hierarchy 
problem requires, at least two extra dimensions are necessarily needed [2j. 

The basic assumption of the extra dimension scenario is that all the fields charged 
under the Standard Model (SM) gauge group are localized on a brane, the familiar 3+1- 
dimensional space-time, embedded in the 4+n-dimensional space, called the bulk |2J. This 
is not only true for the SM non-singlet fermions and scalars, but also for the SM gauge 
bosons, which are attached to the branes e.g. due to a gauge trapping mechanism [3]. In 
string theory the SM particles correspond to open strings whose ends are stuck to the 
branes, while gravitons and other particles that do not carry SM charges correspond to 
closed strings free to propagate in the extra dimensions (see e.g. Refs. [Zl IH] ) - 
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The possible unlocalizability of right-handed neutrinos (vr) in extra-dimension sce- 
narios was emphasized in j^J EH|. This feature may turn out to be of crucial importance 
from the neutrino physics point of view. The Yukawa coupling between the left-handed 
neutrino Vl, residing on the brane, and the right-handed neutrino ur, residing wherever in 
the spacetime, is suppressed by a factor of the order of M^/Mpi due to the scarce overlap- 
ping of the corresponding wave functions. This would offer a new and elegant explanation 
for the lightness of neutrinos. Perhaps the most popular conventional explanation for the 
lightness of neutrinos is given by the seesaw mechanism, which requires the existence of 
a new mass scale around 10 12 GeV or higher. In the extra-dimension scheme such a high 
mass scale is neither needed nor naturally appears but small neutrino masses follow from 
the suppression of the Yukawa couplings due to the large bulk volume. 

The possible existence of bulk neutrinos implies another novel feature. The momenta of 
bulk fields are quantized in the compactified extra dimensions, and each of these quantum 
states corresponds to a massive state, a Kaluza-Klein excitation, at the brane world. 
Hence a bulk neutrino would appear as an infinite tower of massive Dirac neutrinos in 
our four-dimensional spacetime. These new degrees of freedom would hugely enrich the 
neutrino mass spectrum and neutrino oscillation patterns, as discussed in many recent 
papers, see e.g. [TTj. 

In our previous work we investigated a model with one large extra dimension pay- 
ing a particular attention to the effects of the extra dimension in neutrino physics. The 
localization of the SM fermions in the four-dimensional space was achieved by introducing 
a five-dimensional scalar field with position dependent expectation value having a domain- 
wall profile in the extra dimension ^Sl Ej Due to the periodicity of the scalar field in 
the compactified fifth dimension, there necessarily exists, in addition to the ordinary SM 
brane, another brane, a mirror brane, in the bulk. We constructed a concrete mathe- 
matical framework that appropriately described this two-brane scheme. We extracted a 
finite neutrino mass matrix, of which the massive Kaluza-Klein modes were integrated 
out, and studied neutrino masses and mixing from phenomenological point of view. It 
was found out that both the neutrino properties and the size of the extra dimension, R, 
are in this scheme severely constrained by experimental data. Especially, R turned out to 
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be bounded from below due to the upper limit of the electron neutrino mass set by the 
tritium beta decay experiments. The mixing between the ordinary left-handed neutrino 
and the inert left-handed bulk neutrino is fairly small if R ^ 10 -3 mm, but maximal 
mixing becomes possible for a larger R. A large radius, R 10~ 2 mm, seemed also to be 
favored by cosmological arguments. Finally, all the active-sterile neutrino mixing schemes 
considered as possible solutions to the solar and atmospheric neutrino anomalies were also 
found realizable in this model. 

As mentioned before, the hierarchy problem cannot be solved if there is only one extra 
dimension. In this paper we therefore extend our previous analysis to the case of two 
extra dimensions. The two-dimensional extra space is assumed to be compactified on a 
torus. The torus geometry allows for a flat metric, which makes various mathematical 
manipulations quite straightforward. 

We are interested in figuring out in which respect this more realistic six-dimensional 
scheme differs from our previously studied five- dimensional model. In particular, we would 
like to find out how the requirement of sufficiently large extra dimensions - so as to solve 
the hierarchy problem - restricts the parameter space of neutrino physics, that is, mixing 
angles and squared mass differences. 

We will use a similar field-theoretical fermion-trapping mechanism as used in earlier 
studies j2j EE]- We introduce a six-dimensional scalar field that has a vortex-like 
solution, the field vanishing at the center of the vortex The SM fermions, assumed 
to be coupled to the scalar field, are then trapped in the core of the four-dimensional 
vortex, the brane. The periodicity of the scalar field implies that the vortex must have 
two vis-a-vis antivortices in both of the extra toroidal dimensions. Consequently, the 
antivortices require then existence of two vortices, so the vortex configuration consists 
of two vortices and two antivortices. In other words, besides "our" SM brane one has 
another brane where particles have usual chiralities and two mirror branes containing 
mirror particles with chiralities opposite to those of their SM counterparts. In addition to 
gravitons, only right-handed ordinary neutrinos and their opposite-handed mirror partners 
are free to propagate in the higher- dimensional bulk. Note that our way of reasoning here 
is perfectly analogous to the case of one extra dimension ^2] , except that now the number 
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of branes and mirror branes is doubled. 

Note that the vortex structure on a sphere would be simpler with just one vortex- 
antivortex pair. However, the curvature of a sphere is non-zero and therefore no fermion 
zero modes exist on a sphere JB|- As a consequence the bulk (inert) neutrino structure 
and the interplay between the bulk and the brane (active) neutrinos would be much more 
complicated. 

In the following Section we shall introduce the mathematical framework of our six- 
dimensional scheme and study its effective action in the four- dimensional spacetime. The 
neutrino mass matrix is derived and a phenomenological analysis of neutrino physics of 
the model is performed in Section 3. Section 4 is devoted to a summary and conclusions. 

2 Formulation of the model 

In this section we consider a concrete six-dimensional model that exhibits the feature 
of localizing the left-handed neutrinos on branes and the right-handed mirror neutrinos 
on mirror branes while letting their gauge inert opposite-handed companions to live in 
the entire higher- dimensional bulk. The localization of the gauge active neutrinos is 
achieved by the field-theoretical method that was discussed in e.g. Ref. |2] and that 
we applied in our previous analysis of a five- dimensional scenario The two extra 

dimensions are assumed to be compactified on a torus with radii R\ and R 2 - A point in 
the six-dimensional space is denoted as (x^, y\, y 2 ), where x^ stands for ordinary spacetime 
coordinates (// = 0, . . . , 3) and yi( 2 ) ~ 2/1(2) + 27ri?i( 2 ) are periodic coordinates on the torus. 

As mentioned above, the scalar field ($) causing the localization is supposed to form 
two vortices and two antivortices. A vortex solution is given by $ = e ie <f)(r), with 0(0) = 
and <j)(r = large) ~ const., where r is a radial coordinate |17] . For antivortices one simply 
changes the sign of 9. As will be seen soon, neutrinos are trapped inside a vortex, the 
throat of which can be viewed as a four-dimensional brane. 

The neutrino sector of our model consists of six-dimensional chiral spinors, a left- 
handed spinor ^ and a right-handed spinor ^> 2 , both associated with active neutrinos, 
together with a spinor N, which is associated with an inert neutrino and which, as will 
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be seen later, can be chosen right-handed. The fields described by the spinors and ^2 
are confined to branes or mirror branes, whereas N exists everywhere in the bulk. This 
choice is a minimal set-up needed to build a successful model on a torus. However, for 
technical reasons we will keep for a while both chiral components of the spinor fields in 
our considerations. Contrary to the five-dimensional case , two active fields have now 
to be introduced, in order to accomplish the localization of the fields. 
The six-dimensional action of our model is taken to be 



where A = 0,...,5, k is a dimensionful Yukawa coupling, and H corresponds to the 
ordinary Higgs field. The gamma matrices read 



where 7 M 's are the usual four- dimensional gamma matrices, and 7 s = i7°7 1 7 2 7 3 . It is 
easy to show that r A 's fulfil the six-dimensional Clifford algebra {r A ,r B } = 2g AB I 8 , 
with g AB = diag(l, — 1, — 1, — 1, — 1, — 1). Furthermore, by defining the six- dimensional 
chirality operators, P± = ^(I 8 ±T 7 ), where T 7 = — r°r 1 ■ • • T 5 = diag(l 4 , — 14), one is able 
to decompose the six-dimensional spinors as ^1 = + ^i-) T and similarly for \l/ 2 and 
N 1 . The subscript + denotes a right-handed and the subscript — a left-handed spinor 
component in the six-dimensional sense. 

The reader should note that we have not presented the gauge boson contributions, free 
scalar parts, nor the non-neutrino (fermion) parts in the action (JI]), but only the parts 
essential to read out the neutrino couplings. Gauge bosons are, as mentioned earlier in 
this paper, trapped to the vortices, so that, effectively, each vortex has its own SM gauge 
bosons, which couple directly only to fermions on the very same brane. Free scalar parts, 
on the other hand, are the standard ones, and the non-neutrino parts are out of the scope 
of the present paper. 

1 See e.g. Ref. for an earlier discussion on six-dimensional diracology, chiralities in six and four 
dimensions, and other related issues. 



s 
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Let us further point out some peculiar features of our action. The form of the action 
(JIJ follows from the imposition of a global U(l) symmetry, under which the fields 
^2, N and H are taken to have quantum numbers 1, -1, 2, and 1, respectively. 
Especially, as will be seen below, the term ^^i^ 2 takes care of the localization of the 
active neutrinos on the branes. Terms like Q^^i, (i ^ j), H*N^ 2 and 

HN^i are forbidden by the U(l) symmetry. Mass terms and NN as well as the 

bulk neutrino Majorana mass term ~ NN C , allowed by the symmetry, do not survive as 
we will put *i + = = N_ = later on. Combinations $ ( * ) * r i( 2 )^i (2 ) , ^i&j {i ^ j) 
and HN^/2, on the other hand, vanish when the SU(2) structure is taken into account: 
^j's and H are the neutral parts of SU(2) doublets whereas N is a SU(2) singlet. 

We start to analyze our model by studying the localization of the brane fermions. The 
equations of motion of the fields ^2 are 

iT A d A ty 1 + g<5>ty 2 = , 
iT A d A ^2 + g*^i = , (3) 

where a term k*HN has been neglected, since it may be regarded as a small perturbation 
in the lowest order approximation. Using polar coordinates (yi = rcos9,y 2 = rsin^), 
one has 

and thus Eqs. (jHJ) can be written in terms of the chiral components of ^1,2 as 

-i^d^ 1+ = -ie^ e d r ^ 1+ + - 7 5 e^ e 9^ 1+ + ge ie ( p(r)^2- , 

r 

-iYd^ 2 - = te-^ 5e d r ^2- + -l 5 e-^ 9 d e ^2-+9*e- ie (l)(r)^ 1+ , (5) 

r 

-iYd^ 2+ = -ie l ^ e d T ^2+ + -i 5 e^ 5e d e ^2 + + g*e- i0 (j)(r)^i- . 

r 

Here we consider the localization on a vortex (i.e. $ = e te <ft(r)); the case of an antivortex 
localization proceeds along the same lines, with just the chiralities appropriately reversed 
(see below). 
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At this point we return to our original chiral picture by setting = — 0. 
Looking at the first and fourth of equations of (0) 2 , it is now natural to require that 
7 5 ^i_ = and 7 5 ^ 2 + = -*2+- Tlli s means that = |(1 - 7 5 )^i_ = , 

the second sign in the subscript referring to the ordinary four-dimensional chirality of 

(i.e., is left-handed in four-dimensional sense). Similarly \l/ 2+ = \[> 2H , i-e., the 

four-dimensional chirality of ^ is left-handed, while its six-dimensional chirality is "right- 
handed". Furthermore, a solution of the equations of motion apparently exists if de^i- = 
de^2+ = 0. Eqs. (0) now reduce to 

-Z 7 ^*!__ = e"(z0r*i— + <^(r)*2+-) , 

-i-fd^- = e-^-^V + gXr)^). (6) 

Writing = ipi (x)fi(r) and ^ 2 +- — { P2+-(x)f 2 (r), where / 1)2 are real functions, 

and ipi = vip 2+ -, one ends up with the equations 

-A(r)(z 7 ^i-) = eV-(4/iM + ^VM/ 2 W), 
-v~ 1 h{r){i'fd tl <p 1 —) = e^^{-iv-%f 2 {r)+g*<j>{T)h{r)) . (7) 

The l.h.s. of these equations vanish due to the usual four- dimensional Dirac equation 
for a massless field. If we take v = ±z'e* arg< - 9 \ the r.h.s. can be easily shown to have 
solutions behaving like 3 /i, 2 (r) ~ exp(±|g| Jq ds<p(s)), where, in order to confine /i, 2 (r) 
to the vicinity of r = 0, the plus (minus) sign must be chosen for <f>(r) < (> 0). This 
concludes our demonstration of the trapping of the spinor fields ^1 and \I / 2 on a vortex 
by the scalar field $. 

Two comments are worthwhile at this stage. First, it is straightforward to see what 

will happen in the antivortex case. Substituting $ = e~ td (p(r) to Eqs. (JHJ) and (jSJ), one 

sees immediately that now a natural choice is 7 5 \l/i_( 2+ ) = + l I / i-( 2 +), and it is hence the 

right-handed components |_ and \I/ 2 ++ that get localized on an antivortex. This is of 

course what one would have intuitively expected to happen a priori. Second, as can be 

2 The second and third equation would not have entered at all if we had used chiral fields from the 
beginning. 

3 This requires that /i(0) = / 2 (0). 
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deduced from the second and third of equations of ©, the condition = = not 
only makes life easier, but is necessary if only the left-handed neutrinos (ordinary active 
neutrinos) are wanted to get trapped on a vortex and right-handed neutrinos (active 
mirror neutrinos) on an antivortex. As is known, there is so far no indication of the 
existence of mirror neutrinos in our four- dimensional world 
One may now write somewhat symbolically 



». = y/M [ Vl _\ x) ) + y/HV-fk) ( V J_ (X) ) + y/HV-fh) ( V J +(X) 





+ \/%^)(% +W ) , (8) 

where the branes (vortices) are situated at yi(= 0) and y 3 , and the mirror branes (an- 
tivortices) at y 2 an d V&- The form of ^> 2 is due to the requirement that (see above) 
*&i-(r = 0) = ^2+{ r — 0) modulo a phase factor that can be absorbed to $ (or alter- 
natively to g) in the original action of Eq. ((TJ). The issue of the Higgs field is now more 
involved than in the five-dimensional case > and the global behavior of the field is not 
known in general. However, since H can be locally shown to get trapped on a vortex or 
an antivortex [2], we write 



H = sfsi&h^x) + sj5{y-y 2 )h 2 {x) + y/6(y - y 3 )h 3 (x) + ^5{y - y 4 )h 4 (x) , (9) 

where hi can be viewed as the vev of the Higgs field in the zth brane. By substituting 
Eqs. (JHJ) and Q, together with a Kaluza-Klein expansion 

( N + 1 _ 1 V" ( N +( X ) 1 Jkyi/Ri+ily 2 /R 2 fiM 

to Eq. (JTJ), one obtains 

S = J d 4 x {i2^ 1 _„7 M 9 M v9i__ + i2Tp 2 _ + Yd fl (p 2 -+ + i2Tp 3 __^d IJl ip 3 — 
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oo 
k,l=—oo 



+N k L —-—) N k J_ + N k L - + -) 



++ 



\Ri R2 ) \ -Ri R2 , 

+u{h\N k l + ^ + h* 2 z k 2 l N+_V2- + (11) 
+h* 3 z kl N+ + ip 3 — + hlzfN k l_ip^ + + h.c.)] } , 



where 



u = , (12) 

k is taken to be real, and z kl = e - '^'"'^^ (t^-, t^-) indicating the position of the j'th brane 
in the extra dimensions (with $1 = oui = 0). Our analysis complies with a "minimalistic" 
view in the sense that we have set iV_ = 0, i.e., N is taken to be a right-handed chiral 
field. We remind that this choice made the Dirac and Majorana bulk neutrino mass terms 
to disappear in the original action ([TJ. 

In order to integrate out the massive Kaluza-Klein modes, we proceed along the same 
lines as we followed in our previous study [T2]. This time, however, the ensuing equa- 
tions are quite lengthy and intransparent and they are mostly omitted in the following 
presentation. The first ten k, /-dependent terms of the action (fTTj) disappear due to the 
equations of motion of the Kaluza-Klein excitations N kl , for example, 

N k + l + = - f A + JjQ (*7^+- + V+ + «^4W+) , (13) 

and similarly for N kl _ (here obviously (k,l) 7^ (0,0)). Substituting these twice to the 
remaining four h.c. -terms of the action, neglecting higher order derivatives of the bulk 
fields, and regarding h^s as constants, one obtains 

£(M)tKo,o) = -i^RxiE^Lpi) + h.c.) + iu 2 pRi^^ifi) - h.c.) 

+m 2 E?[|/ il | 2 S (0,0,p)^ 1 __ 7 ^ 1 __ + |/i 2 | 2 S (0,0,p)<p 2 _ +7 ^ 2 _ + 
+ |/i3| 2 s(0,0,p)^3„_7 M <9^3-- + |/i 4 | 2 .s(0,0,p)^4_ + 7 At <9 At v?4-+ 

+ (/tl/l3s(tf 3 l, W31, p)^!_„7 M ^^3— (14) 

+h 2 hls(§ 42 , ^42, p)Tp 2 _ + ^d^ A - + + h.c.)] , 
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where i?y = $j — cjjj = cjj — Wj, 



S* (w) (0,W,p) = k 2 + n 2]2 e 

(fc,0^(o,o) * 

Hu,(tf)(<Pi) = (+) h 2 h\s ^012, Wi2,p)^ 2 -+V 9 l 

(+) /l 4 ^l^(o;)(^14,^14,P)^4- + V 2 l — 
+ /l 3 /l2S^)(^23, ^23, p)V 3 — ( P2-+ 
(+) /i4^3 S «?H(^34,^34,P)^4-+V 9 3— , 

and p = R1/R2. Combining then Eqs. (fTljl and (fTT|) (with = I — 0), we end up with 
the effective Lagrangian 

C e ff = iTp 1 ^d^(p 1 + iTp 2 Yd^(p2 + i^ 3 ^d^(p 3 + ilp 4 Y9 l ,(p4 + iN + 'y fl d f ,N + 
+u(h* 1 n 1 N ++ (p 1 + h* 2 n 2 N + _ip 2 + h%n 3 N ++ Lp 3 + h* 4 n 4 N + _(p 4: + h.c.) 
-u 2 Ri(E&(ni<pi) + h.c.) + iv 2 pRiCEuiniipi) - h.c.) (16) 
+i(cffp 1 'fd lx (p$ + bip^d^ipi + h.c.) , 

where we have suppressed the chirality indices of <pj's and the Kaluza-Klein indices of 
iV's, rescaled <pj fields with rij = (2 + u 2 R 2 \hi\ 2 s(0, 0, p))~ 1//2 , and we have denoted 

o = u 2 R\h 1 hlnin 3 s(d 3 i,uj 3 i, p) , 

b = u 2 Rlh 2 hln 2 n 4 s(i!} A 2^A2, p) ■ (17) 

Note that s(0, 0, p) is formally a divergent function. Therefore physical cut-offs k < RiM*, 
I < R2M* have been introduced in the sums. This leads to the expression 

, Mpi 7T 2 / 1 \ 2/1 1 

p S (0,0,p).2,ln— + . + T (p + -) --(-+- 
where M* = (M 2 J R^) 1 ^ . 
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Finally, by denning 
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(19) 

in order to get rid of the kinetic cross terms, and rescaling £j's with Wi(3) = (1(+) |a|) _1 ^ 2 
and n 2 (4) = (1(+) |6|) -1 / 2 , Eq. (|TE)l may be cast into the form 

-(Al 2 ^ + BU, + CU, + D£ 4 £ 3 + h.c.) (20) 
-(£iV ++ 6 + + GN++& + HN + M + h.c.) , 



where 



£ 
T 

Q 
H 



'72 
un 2 

un 3 

wn 4 
'72 



(/i*nie iCa - hln^e~ iC ' a ) , 
(^n 2 e iC " - /i^e - * 6 ) , 
(/i^ie iCa + h* i n ? ,e~ iC - a ) , 
(^n 2 e^ + hln A e~ i(b ) . 



(21) 



The explicit expressions of the quantities A, B, C and T> are irrelevant for what follows. 
One should notice that if the phases of /ij's are absorbed to the <pi fields in Eq. (ITTj) . one 
has, s(#, p) being real, Ca = Cfe = 0, and thus also 8, J 7 , Q and H can all be made real. 

3 Neutrino phenomenology of the model 

We now move to consider the phenomenological aspects of our scheme, in particular 
the effects of the multibrane structure on neutrino physics. Apart from gravitons, only 
sterile neutrinos can mediate contacts between branes and mirror branes through the 
bulk, thereby affecting neutrino phenomenology at our home brane. 
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From the effective Lagrangian £ e ff m Eq. (|2U|) one infers the following mass Lagrangian 
for neutrinos: 



C M = -{i 2 UN ++ )M 
where the matrix M is given by 

M = 




+ h.c. 



(22) 




(23) 



Even though M can in principle be diagonalized in its general form, this task becomes 
more straightforward if the branes are assumed to be distributed symmetrically on the 
torus 4 . We have made this assumption in the following. Then, as can be deduced from 
Eqs. (JUJ), H# = 2^ = and consequently A = B = C = V = 0. As will be seen later, 
even this simplified scenario leads to diverse neutrino phenomenology. Defining 



O a 



( 
\ 



cos a sin a 
- sin a cos a 
1 



with 



we obtain 



cos a 



V£ 2 + Q 2 



cos j3 



( cos/? sin/5 
— sin (3 cos (3 
1 



H 



VF 2 + H 2 



O^MO a = diag(mi,m 2 ,m 3 ) 



(24) 

(25) 

(26) 
The corre- 

(27) 



The neutrino sector thus consists of one massless and two massive Dirac neutrinos. We 
emphasize that the existence of one massless eigenstate is not an intrinsic feature of our 



where the eigenvalues are mi 
sponding mass eigenstates are 



0, m 2 = V£ 2 + Q 2 and m 3 = ^ T 2 + H 2 



(Xi \ 




X2 




V X3 j 






4 For example -#3 = ^3 = t?2 = ^4 = tt, ^2 = $4 = 0. Then s(i?3i, W31, p) = s(i?4 2 , 0J42, p) — s(tt,tt,p), 
which is easily calculable numerically. 
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scenario but results from our choice of symmetrically located branes. Another charac- 
teristic feature of the model is the imposed right-handed chirality of the bulk field N. 
As mentioned before, this choice automatically wipes out the six-dimensional Dirac and 
Majorana bulk neutrino mass terms, thus removing a usual need to rely on ad hoc argu- 
ments. 

Let us remind the reader that in Eq. (}2"T)) £1 and £3 are related to left-handed active 
neutrinos (through e.g. Eq. (|19jl ). £2 and £4 to right-handed active mirror neutrinos and 
the N fields represent sterile bulk neutrinos not confined to any brane. Using Eqs. (|19p. 
(|24[) and (J2"7)l . taking into account the field rescalings, and normalizing correctly, the 
original left-handed neutrino of "our" brane is then given by 



As ipi is the only active neutrino living on our brane, only the mixing angle 9 can be 
experimentally probed. The angle /3 is related to the right-handed fields <pi and ^4, and 
is therefore a measurable quantity only on the mirror branes. 

It should be pointed out that the physical setting of the neutrino mixing differs in 
one noteworthy respect from the results of our previous five-dimensional study. In 
the ordinary left-handed neutrino was found to be mixed with an inert left-handed bulk 
neutrino. Now, as can be seen easily by repeating the calculations similar to those leading 
to Eq. (J2EJ), the other linear combination (not necessarily orthogonal to Eq. (J2BJ)) of Xil 
and X2L is ips- fi is hence mixed with a left-handed active neutrino living on another 
brane. 

Similarly, as seen from Eq. (|28j) . the massless neutrino xi an d the massive state X2 
appear as mixed in weak interactions. There exist several empirical constraints on their 
mixing angle 9 and the mass 777.2, coming from laboratory experiments, astrophysical 
observations and cosmological considerations. In Fig. 1 we compile (as a rough step- 
function approximation sufficient for our purposes) the existing laboratory bounds on 



tp x = cos 9xil + sin 9x2L , 




where 



cos 9 



rii c °s a — n 3 sin a 



(29) 
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the mixing angle 9 between the electron neutrino and a heavier neutrino in various mass 
ranges of the mass 7712, see Refs. |2T]. These bounds follow mainly from the lack of 
extra kinematical thresholds in the spectra of various particle decays (e.g. beta decay or 
7r + — > e + v e ) or the nonobservation of the heavier neutrino decays like v — > u e + e~ + e + . 
One should take into account also the upper limit of the electron neutrino mass from the 
tritium beta decay experiment 



m v „ <> 2.2 eV. 



(30) 



This bound should be fulfilled by the mass eigenvalue rri2, whenever the electron neutrino 
is predominantly the heavier mass eigenstate \2- An important constraint at low 7712 
values is given by the oscillation limit from the Bugey disappearance experiment j2S]- 
The upper bound for sin 2 29 varies in the range 0.02 to 0.3 for m 2 , in the range 0.01 eV 2 
to 100 eV 2 . 



sin 9 




Figure 1: Laboratory bounds on the mixing angle between the electron neutrino and a 
sterile neutrino from laboratory measurements as a function of the mass of the heavier 
mass eigenstate (the lighter eigenstate is massless). 



The active-sterile neutrino mixing is constrained also by cosmology. If the mixing is too 
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large, neutrino oscillations, acting as an effective interaction, would bring sterile neutrinos 
into equilibrium with the light SM particles before neutrino decoupling. The resulting 
excess in energy density would endanger the standard scheme of the nucleosynthesis of 
light elements [2I|. This leads to the following bound for v e <-> v s mixing 25 : 



(This bound is avoided if there was a suitable net lepton number in the early universe 



We will now confront our scheme with these constraints. We will search the allowed 
regions of the parameter space of the model numerically using a Monte Carlo analysis. 
We vary three unknown parameters: the radii of the two extra dimensions R\ and R 2 , and 
the Higgs vacuum expectation value /13. The extra dimension radii are varied from the 
Planck scale to a millimeter scale and the Higgs vev between 10~ 2 MeV and 10 12 MeV. 
Both i?j's and h 3 are randomized so that their logarithms are evenly distributed. For h\ 
we use the value 174 GeV, and, as can be seen from Sees. 2 and 3, h 2 and /14 do not affect 
neutrino physics on our brane. 

In addition to these parameters the action depends on the dimensionful Yukawa cou- 
pling k. The natural scale for it can be thought to be set by the higher-dimensional 
Planck scale M*, and therefore we write (with Mp t = M^RiR 2 ) 



where k' is a dimensionless coupling constant. The value of k' is not really known, but 
a plausible choice would be a number relatively close to unity. We have performed our 
analysis for the values k? — 1 and k' = 0.001. 

In Figs. 2 we present scatter plots in the (R, m2)-plane (with R = \fR~iR~2) for the 
two values of the Yukawa coupling (/«' = 0.001 and 1). These figures consist solely of 
"raw data", i.e. neither the cosmological constraint, Eq. (|3*T|) . the limits of Fig. 1, nor 
the Bugey constraint are taken into account. The allowed parameter space is seen to be 
quite severely constrained by the mathematical structure of the model. Nevertheless, the 
characteristic extra-dimensional size may vary from the millimeter scale some 30 orders 



\5m 2 \ sin 2 29 < 5 x 10~ 8 eV 2 for \5m 2 \ < 4 eV 2 , 
sin 2 26 < 10~ 8 for \5m 2 \ > 4 eV 2 . 



(31) 




(32) 
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of magnitude all the way down to the Planck scale, and the higher mass eigenvalue m 2 
ranges from 10 -9 eV to 10 10 MeV, depending on the coupling k'. The general trend is 
that small (large) extra-dimensional size corresponds to large (small) mass. Note that 
the value of k' affects the vertical position of the two diagonal borderlines, whereas the 
wedge-like region below them remains immutable. 




Figure 2: Scattering plot of geometric average of the radii, R, vs. possible neutrino masses 
m 2 without any experimental constraints, a) k' = 0.001, b) k' = 1.0. 

Figs. 3 show how the situation changes when the experimental constraints are taken 
into account. Interestingly, the extreme values on both axes are still allowed and, vaguely 
speaking, only dots in between are cut out. Numerical investigations suggest that the 
upper (lower) diagonal borderline reflects the chosen range for the Higgs value h 3 , higher 
(lower) values enabling mass m 2 to be still larger (smaller). For the lower diagonal sin 9 
1. In the case of k' — 1 the allowed region is split up into two separate parts, the feature 
that will be seen in all the subsequent figures. 
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Figure 3: Scattering plot of geometric average of the radii, R, vs. possible neutrino masses 
m2 with all observational constraints, a) k' = 0.001, b) k! = 1.0. 

Figs. 4 display the behavior of the mixing angle 9 as a function of R. For very 
small extra-dimensional sizes the mixing is tiny, but above R ~ 10~ 16 mm maximal 
mixing, 9 = ir /4, becomes possible. The peculiarity that 9 can never be exactly zero is 
a consequence of the mathematical structure of the model, as may be deduced from Eq. 
(J2HJ) . However, for larger values of /13 the allowed region would extend to smaller angles. 
The limit just below sin 9 = 10 -4 , on the other hand, is due to the cosmological constraint 
(f3*T]) . Note that the obvious change in the scatter plot for R <: 10 -10 mm (V = 0.001) or 
R <: 10~ 4 mm (k' = 1), as well as the tilted rightmost borderline penetrating both of the 
separate parts in the case of k' — 1, have analogous counterparts in Figs. 3. 

Finally, in Figs. 5 is shown the allowed parameter space in the (sin 2 29, <5m 2 )-plane, 
where 5m 2 = m\ as m\ = 0. We have presented in these figures separately those points 
which correspond to large values of the size R of the extra dimension, 0.01 mm < R < 1 
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Figure 4: Scattering plot of geometric average of the radii, R, vs. mixing angle sin 9 with 
observational constraints, a) k' = 0.001, b) k' = 1.0. 



mm. With R in this range, the hierarchy problem can be solved in this model and also the 
existence of the extra dimensions is testable in the new short-range gravity experiments 
[HI EZj • The empty areas in the large sin 2 29 large 5m 2 region are forbidden by the 
cosmological constraint The other constraints do not have substantial effects to this 
plot. As mentioned before, the borderline at small 9 (now at sin 2 29 ~ 10 -13 ) corresponds 
to the upper end of the Higgs vev h 3 range covered by our Monte Carlo data set, and it 
would move towards smaller values if the maximum value of the parameter is increased. 

As can be seen from Figs. 5, our model allows for a vast diversity in the possible 
(5m 2 , ^-combinations. With the chosen ranges of the free parameters of the model, the 
squared mass difference of the active-sterile neutrino pair can have values from about 
10~ n eV 2 (V = 1) or 10 -17 eV 2 (V = 0.001) up to 10 25 eV 2 or more. However, squared 
mass differences Sm 2 = m\ ^ 1 eV 2 are possible only for very small values of the mixing 
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angle, 9 ^ 10^ 4 . If the extra dimension size is in the theoretically and phenomenologically 
interesting region of R <; 0.01 mm (the hits denoted by circles in Figs. 5), 5m 2 is forced 
to much smaller values. In particular, maximal active-sterile mixing is possible in this 
case only for quite small values of the squared mass difference: 5m 2 ~ 10 -11 — 10~ 9 eV 2 
for k' — 1 and 5m 2 ~ 10~ 17 — 10~ 15 eV 2 for k' = 0.001. On the other hand, the largest 
squared mass difference consistent with the large extra dimension size is of the order of 
5m 2 ~ 10~ 6 eV 2 . 

Let us still mention that the leftmost points in the small squared mass difference region, 
5m 2 < 1 eV 2 , correspond to 9 ir/2, as may be seen particularly well by comparing Figs. 
4b and 5b. In this case the predominantly active state, the ordinary neutrino, is the 
heavier one of the two mass states. 
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Figure 5: Scattering plot of sin 2 (2#) vs. squared neutrino mass difference Sm 2 with ex- 
perimental and observational constraints taken into account. Circles indicate the points, 
where 0.01 mm < R < 1 mm. a) k' = 0.001, b) k' = 1.0. 
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4 Summary and conclusions 



We have investigated a six-dimensional model with two extra dimensions compactified 
on a torus. Particles with Standard Model charges are trapped on four four- dimensional 
branes (vortices) by a scalar field with an appropriate periodic position dependence on 
the torus. On two of the branes, one of which is the world we live in, particles have 
the familiar V — A weak currents, while on the other two, called mirror branes, weak 
currents are of the V + A form. The structure of the model is mathematically quite 
strictly determined, but there are some free parameters in the theory such as the size of 
the extra dimensions, characterized by the radii of the torus R\ and R2, and the vev (h 3 ) 
of the counterparts of the SM Higgs field on the other branes. Also the strength k' of the 
Yukawa coupling between the Higgs field, brane neutrinos and bulk neutrinos, as well as 
the locations of the vortices on the torus, are not a priori restricted. 

Right-handed neutrinos and left-handed mirror neutrinos, which lack Standard Model 
interactions, can propagate in the extra two-dimensional bulk and thereby mediate mix- 
ing between active neutrinos of the branes and mirror branes. This affects neutrino phe- 
nomenology on our brane. We have confronted the model with existing data on neutrino 
masses and active-sterile neutrino mixing from laboratory measurements and cosmology. 
Our main results can be seen in Figs. 5, where we have plotted the (sin 2 29, 5m 2 )-parameter 
space covered by our model when the free parameters of the model are allowed to vary 
within wide ranges of values. In order the model to solve the hierarchy problem, the size 
of the extra dimensions should be large enough. The corresponding points in Figs. 5 are 
denoted by circles, and as one can see, the squared mass difference is in that case forced 
into a quite narrow range of values, in particular when the mixing is maximal. 

We have restricted our phenomenological analysis to the electron neutrino, but con- 
clusions would be very similar for the muon and tau neutrinos, since the most important 
constraint, the cosmological bound for the active-sterile mixing, is quite insensitive to 
neutrino flavor. 

Let us finally emphasize that the model we have considered is not the most general one 
but we have made a number of simplifying, albeit conceivable, assumptions. Nevertheless, 
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we believe that our analysis has revealed the essential features of the scenario. 
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